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ABSTRACT 

We use numerical simulations to study the relation between the velocity of the Local 
Group (LG) and its gravitational acceleration. This relation serves as a test for the 
kinematic origin of the CMB dipole and as a method for estimating f3. We calculate the 
misalignment angle between the two vectors and compare it to the observed value for the 
PSCz survey. The latter value is near the upper limit of the 95 % confidence interval 
for the angle; therefore, the nonlinear effects are unlikely to be responsible for the whole 
observed misalignment. We also study the relation between the amplitudes of the LG 
velocity and gravity vectors. In an Q = 1 Universe, the smoothed gravity of the LG turns 
out to be a biased low estimator of the LG (unsmoothed) velocity. In an Q — 0.3 Universe, 
the estimator is biased high. The discussed biases are, however, only a few per cent, thus 
the linear theory works to good accuracy. The gravity-based estimator of the LG velocity 
has also a scatter that limits the precision of the estimate of /3 in the LG velocity-gravity 
comparisons. The random error of j3 due to nonlinear effects amounts to several per cent. 



1 Introduction 

The dipole anisotropy of the Cosmic Microwave Background (CMB) temper- 
ature is widely believed to reflect, via the Doppler shift, the motion of the Lo- 
cal Group (LG) with respect to the CMB rest frame. When transformed to 
the barycenter of the LG, this motion is towards (I, b) = (276° ±3°, 30° ±2°), 
and of amplitude ut_g = 627 ±22 km-s" 1 , as inferred from the 4- year COBE 
data (Lineweaver et al. 1996). Alternative models which assume that the 
dipole is due to a metric fluctuation (e.g., Paczyhski & Piran 1990) have 
problems with explaining its observed achromaticity and the relative small- 
ness of the CMB quadrupole. 

An additional argument in favour of the kinematic interpretation of the 
CMB dipole is its remarkable alignment with the LG acceleration, inferred 
from galaxy distribution. The acceleration on the LG (i.e., the galaxy 



1 



2 



dipole), inferred from the IRAS 1.2 Jy survey points ~ 25° away from the 
CMB dipole (Strauss et al. 1992; hereafter S92). The recently completed 
IRAS PSCz survey allowed to make further progress on this topic. Schmoldt 
et al. (1999; hereafter S99) found the PSCz galaxy dipole to be within 15° of 
the CMB dipole. Rowan- Robinson et al. (1999), performing a similar anal- 
ysis of the galaxy dipole but out to a larger distance 300 h^ 1 Mpc, obtained 
the misalignment angle as small as 13°. 

In the linear theory, (e.g., Peebles 1980) the peculiar velocity of any 
galaxy, v, is directly proportional to its gravitational acceleration, g: v = 
(2f /3HQ)g . Here, ft denotes the cosmic matter density parameter, H is 
the Hubble constant and 

/(n,A)~n°- 8 + A^i + ^ (i) 

(Lahav et al. 1991; note a very weak dependence on the cosmological con- 
stant). 

The acceleration is caused by the gravitational pull of the surrounding 
matter, 

g(r) =G Pb J dV<5(r')|^3- (2) 

Here, G is the gravitational constant, p b is the background density, and 5 
is the mass density fluctuation field. If we define the scaled gravity, g = 
(2//3-ffT2)g, then in the linear theory 

v = g. (3) 

The vector g can be measured from redshift surveys which give an estimate 
of the three-dimensional galaxy density field, 8 g . In the simplest biasing 
model the galaxy and mass distributions are linearly related, 5 g = bS, where 
b is a linear biasing factor. The scaled gravity is then 

where (3 = f jb and we express distances in units of km • s _1 . 

The amplitude of the scaled gravity depends on (3. Therefore, a com- 
parison between the LG scaled gravity and the CMB dipole can serve not 
only as a test for the kinematic origin of the latter but also as a measure of 
the parameter (3. Combined with other constraints on bias, it may yield an 
estimate of itself. 
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However, the estimate (Q) of the LG velocity from a particular redshift 
survey will in general differ from its true velocity, for a number of reasons. 
We enumerate them here following S99: 

• The finite volume of the survey - the galaxy dipole may be affected 
by contributions from outside; 

• Unsurveyed regions within the volume; 

• The shot noise - the galaxy density field is sampled discretely, at 
galaxy's positions; 

• Redshift-space distortions - in redshift surveys the radial coordinate 
is not the true distance to a galaxy but rather a sum of the distance 
and the radial component of its a priori unknown peculiar velocity; 

• Nonlinear effects - the scaled acceleration is equal to the velocity only 
in the linear regime, while redshift surveys unveil nonlinear structure 
of the galaxy distribution. 

In the proper process of the LG gravity-velocity comparison, all these effects 
should be accounted for. In the present paper, we will concentrate on the 
nonlinear effects (hereafter NE). In general, the NE can be due to both 
nonlinear gravity and nonlinear biasing. Here, we will only consider the 
effects of nonlinear gravity, and we will use the term 'NE' in this, more 
narrow, meaning. 

NE modify the linear relation (||) in a number of ways. Non-local nature 
of gravity, somewhat hidden at linear order, manifests itself at higher orders. 
The relation between the velocity and the scaled gravitational acceleration 
becomes not only non-linear but also non-local, so at a given point it has 
a scatter. The NE may also spoil the alignment between the two vectors. 
Thus, the NE may influence not only the accuracy of estimating [5 from 
Local Group dipole comparisons, but also its precision. In this paper we 
address both issues. 

A common method of constraining cosmological parameters by the LG 
gravity-velocity comparison is to apply a maximum-likelihood analysis. In 
a given cosmological model, one maximizes the likelihood of measuring the 
scaled gravity of the LG given the true (CMB-inferred) value of the LG 
velocity (Juszkiewicz, Vittorio & Wyse 1990, Lahav, Kaiser & Hoffman 1990, 
S92, S99). In such an analysis, a proper object describing the nonlinear 
effects is the decoherence function, i.e. the cross-correlation coefficient of 
the Fourier modes of the gravity and velocity fields (S92). 
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Here, we will study the NE in real space. In particular, we will study the 
evolution of the misalignment angle, and the relation between the amplitudes 
of the gravity and velocity vectors. This approach is more appropriate in 
the case of the 'numerical analysis' of S99, where they simply equated the 
LG gravity, inferred from the IRAS PSCz catalog, to the LG velocity. To 
account for the NE (and other effects mentioned earlier), S99 used mock 
catalogs to compute the ratio between the reconstructed gravity dipole at 
the observer's position and its true N-body velocity. Then they used the 
average of the ratios calculated from mock catalogs as a multiplicative factor 
which relates the reconstructed gravity dipole to the real LG velocity. The 
mock catalogs were constructed for different cosmologies. In our opinion, 
if the factor varies systematically from model to model, averaging it makes 
no sense. Rather, the results should be expressed explicitly in a model- 
dependent way. (This is the case of the likelihood analysis of S99.) One of 
the goals of the present paper is to clarify this issue. 

The NE in real space were studied by means of N-body simulations by 
Davis, Strauss & Yahil (1991); however, only in a standard CDM model. 
Here we investigate the effects of varying on the NE, also using numerical 
simulations. Instead of using a N-body scheme, we model cold dark matter 
as a pressureless cosmic fluid (see Peebles 1987). The outline of this paper 
is as follows: In Section || we present our numerical model of the large 
scale structure evolution and how we select Local Group candidates in the 
data. Next we discuss the differences between the Local Group velocity and 
acceleration: in Section || we discuss the angle between the two vectors, in 
Section || - their amplitudes. We summarize our results in Section [|. 

2 The simulations 

2.1 The Numerical Model 

We have performed our simulations using an Eulerian, uniform-grid based 
code - CPPA (Cosmological Pressureless Parabolic Advection, see Kudlicki, 
Plewa & Rozyczka 1996, Kudlicki et al. 2000a, Kudlicki et al. 2000c). The 
main features of CPPA are parabolic density and velocity profiles, variable 
timestep, periodic boundary conditions and a flux interchange procedure, 
implemented as an approximation to the solution of the Boltzmann equation. 
We have studied two cosmological models with £1 = 1 and f2 = 0.3, assuming 
Gaussian initial conditions. For better statistics we performed 4 realizations 
of each of the models, varying random phases of the initial density field. The 
grid we used was 64 3 , and the comoving size of the simulation box was 100 
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h Mpc on a side. 

The linear velocity depends on the cosmological constant very weakly 
(see eq. |]]); this holds also for higher orders (see Appendix B.3 of Scoccimarro 
et al. 1998 and Nusser & Colberg 1998), so we were free to assume A = 
in all our models. To make the simulated gravity of the LG as close as 
possible to that inferred from the IRAS PSCz survey, for the mass power 
spectrum we adopted the power spectrum of the PSCz galaxies (Sutherland 
et al. 1999): 

B± (5) 

{1 + [ok + (6A;) 3 / 2 + {ckff) 2 ^} K ' 

6.4/r/i -1 Mpc, b = 3.0/r/i _1 Mpc, 
1.7/r /i^Mpc , i/ = 1.13 

with r = 0.2 as best fitted value. To normalize the power spectrum we used 
the observed local abundance of galaxy clusters. The present value of ag, 
(78,0, is a function of f2 and for the case of A = 0, considered here, it is given 
by the relation (Eke, Cole & Frenk 1996) 

a 8 ,o = (0.52 ± 0.04)ir°' 46+0 ' 10Q . (6) 

This relation changes only slightly with the shape of the power spectrum. 
It is also very similar for the case A = 1 — f2. For f2 = 1, as t o — 0.52, while 
for n = 0.3, ct 8i0 ~ 0.87. 

2.2 Selection of LG candidates 

To study the relation between the velocity and the gravity of the LG, from 
the simulation grid we selected the cells which have properties resembling 
those of the LG. We chose these 'LG candidates' following the criteria pre- 
sented in Gorski et al. (1989), and used by Davis, Strauss & Yahil (1991): 
the peculiar velocity of the candidate must be |v| = 600 ± 10 km • s , the 
candidate must be in a region in which the fractional overdensity, averaged 
in a radius R = 500 km • s -1 , S R , is in the range -0.2 < 5 R < 1.0, and 
the cosmic flow within this volume is quiet: |v — v R | < 0.3|v|, where v R 
is the mean velocity within the averaging sphere. In general, the symbol 
x R denotes the quantity x averaged with a top-hat filter of radius R. The 
chosen value 600 km • s _1 is not exactly equal to the present estimate of the 
LG velocity, given in Section ||. However, for the purposes of the present 
paper the difference between the two is not significant, while the value 600 
km • s _1 facilitates a visual analysis of the results. 
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3 The misalignment angle 
3.1 Analytic predictions 

As stated in Section |l], the NE spoil the alignment between the velocity and 
gravity vectors. In this Section we study the evolution of their misalignment 
angle, a, given by 

cos a = . (7) 

vg 

In the linear regime, v is proportional to g and the misalignment angle is 
exactly zero. The linear velocity depends on f2 (and only weakly on A) solely 
via the multiplicative factor / (eq. [l]). This is also approximately true for 
the nonlinear velocity (Appendix B.3 of Scoccimarro et al. 1998 and Nusser 
& Colberg 1998). Moreover, the gravity field scales linearly with Q (see 
eq. [2]). As a result, in equation (|7|) the factors / and Q cancel out and we 
expect the misalignment angle to be practically independent of Q. 

If the fields are smoothed on mildly nonlinear scales, the perturbation 
theory can be applied to predict qualitatively the evolution of a. Note first 
that the angle between the velocity and the acceleration equals to the angle 
between the velocity and the scaled acceleration (eq. ^). Expanding the 
velocity and the scaled gravity in perturbative series, v = Vi + V2 + . . ., 
g = gi + g2 + . . ., and calculating the expectation value of the angle to the 
leading order we obtain 

/ 2\ / (V2 ~ g2) 2 ~ [Vl ' (V2 ~ g2)] 2 

{a ) = ( - 2 

\ n 

Here, vi = vi/vx, and the symbol (• • •) denotes ensemble averaging. 
Perturbative solutions for V2 and g2 (e.g., Goroff et al. 1986) yield the 
following approximate scaling relations: v 2 ~ c/2 ~ ui^i- Since to the leading 
order 81 = S and typically 5 ~ a, we have 

(a 2 ) 1/2 (xa, (9) 

where a = (<5 2 ) 1//2 is the r.m.s. fluctuation of the mass density field. As 
described earlier, we expect the coefficient of proportionality in the above 
relation to be insensitive to O. On the other hand, we expect it to depend on 
the relative amount of power on small (nonlinear) scales. In equation (||), a is 
expressed in terms of only first and second-order perturbative contributions 
to the velocity and gravity fields. Despite this, an analytic calculation of 
the misalignment angle is impossible, because it involves averaging of a 
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ratio of two non-Gaussian variables. In short, the resulting series cannot be 
truncated. 

Thus far, our analysis concerned the evolution of the misalignment angle 
of the smoothed velocity and gravity fields. In practice, the reconstruction 
of the LG gravitational acceleration from the redshift-space galaxy field does 
involve smoothing (to mitigate strong nonlinear effects, to reduce shot-noise 
and distance uncertainty, etc.). The LG velocity, however, is inferred directly 
from the CMB dipole, so it is not a subject to any smoothing, except the one 
needed to reflect the finite size of the LG. As that smoothing, S92 adopt a 
top-hat of radius 100 km-s -1 (1 Mpc). Thus, while the LG acceleration 
is smoothed, the LG velocity is (almost) unsmoothed. Though in the linear 
regime = g^, if v 7^ v K , then v 7^ g R . In other words, if some part of the 
central velocity is due to matter distribution within the smoothing radius 
(what is quite natural to expect) we cannot expect unsmoothed velocity 
to equal smoothed acceleration even in the linear regime (see also Berlind, 
Narayanan & Weinberg 2000). 

3.2 Numerical results 

We have calculated the misalignment angle from our simulations. First 
we smoothed both velocity and gravity fields with a top-hat of radius 500 
km ■ s~ 1 . This radius is equal to the minimum radius of smoothing used 
for the IRAS gravitational field calculations. As a characteristic angle we 
adopted the quantity (a 2 ) 1 / 2 , which, for simplictity, we will denote a. 

Figure Q shows the time-evolution of a for all points of the simula- 
tion as a function of the r.m.s. value of mass fluctuations in spheres of 
canonical radius 800 km • s _1 , erg. The angle depends linearly on as up to 
ag — 0.5. Moreover, it is practically independent of ft: the systematic dif- 
ference between the values for models with the same phases of the initial 
inhomogeneities and different f2 is not greater than the scatter for models 
with the same O and different phases. 

Our second step is to apply different smoothings to the final gravity and 
velocity vectors, to mimic better real situation. The gravitational accelera- 
tion is smoothed with a filter described earlier, but the velocity should now 
be smoothed with a top-hat of radius 100 km • s^ 1 (to reflect the finite size 
of the LG). Since the cell volume in our simulations corresponds to that of 
a sphere with radius of 97 km ■ s , the unsmoothed velocity field is used for 
simplicity. 

Figure [2] shows distributions of a for different models with the above 
normalizations. The mean values of the characteristic angle are 7.0° for 




Figure 1: The characteristic angle between the smoothed velocity and 
smoothed acceleration of all points, for models with $7 = 1 (solid lines) 
and O = 0.3 (dotted lines), as a function of the density dispersion in spheres 
of comoving radius 800 km • s _1 . 



Q = 1 and 7.7° for Q = 0.3. On the confidence level of 95 %, a is smaller 
than respectively 12.6° and 13.9° in the model with f2 = 1 and Q = 0.3. 
Thus, the measured misalignment angle between the LG velocity and PSCz 
galaxy dipole (13-15°) is rather unlikely to be caused entirely by nonlinear 
effects. Other effects must play a role here as well, like, e.g., incomplete sky 
coverage and shot noise. 

The values of the angle we obtained are smaller than the value obtained 
by Davis et al. (1991) from N-body simulations, 10.2°. However, we use the 
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Figure 2: The simulated PDF of the misalignment angle between velocity 
and smoothed gravity for LG candidates. Upper panel: £1 = 1.0; lower: 
f2 = 0.3. Shaded area represents observational constraints from the PSCz 
catalogue. 



spectrum of the PSCz galaxies, while the power spectrum they used was 
a standard CDM. This model has more power on small (nonlinear) scales, 
so we expect a larger effect in this case. We additionally calculated the 
angle in a standard CDM model (the power spectrum given by eq. [| with 
r = 0.5). The mean angle indeed increased from 7.0° to 8.1°. The rest 
of the difference we attribute to differencies between the numerical codes 
(hydro vs. N-body). 
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4 The relation for amplitudes 

In the previous Section we have shown that the misalignment angle between 
the vectors of velocity and gravitational acceleration is small. Therefore, 
most information about the relation between the two vectors is contained 
in the relation between their amplitudes. This is also the crucial point in 
determining Q from the LG dipole. Figure || shows this relation for the LG 
candidates in models with Q = 1 and O = 0.3. Here, both gravity and 
velocity fields are smoothed with a top-hat filter of radius 500 km • s _1 . 

In the top panel of Figure p| the points agree with the linear prediction 
(solid line) quite well. The reason of this is twofold. First, the model 
with f2 = 1 has a low normalization, so the velocity and gravity fields are 
only weakly nonlinear. Second, the velocity 600 km • s^ 1 is rather typical 
for an SI = 1 universe, while the velocity-gravity relation deviates more 
significantly from the linear prediction only in the high- velocity tail (Kudlicki 
et al. 2000b). Specifically, the r.m.s. value of the velocity of all points in the 
simulation is 458 km • s" 1 . 

Contrary to the previous case, in the lower panel of Figure ^ the points lie 
farther off the line. This is not surprising, because the normalization of the 
model with $7 = 0.3 is higher, so the fields are more nonlinear. Moreover, 
the LG velocity of 600 km • s _1 is less typical for an = 0.3 Universe. 
Higher normalization of the model compensates for this effect only partly. 
The r.m.s. value of the velocity of all points in this model is 340 km • s _1 . 

The above discussion does not, however, imply that the gravitational 
acceleration of the LG is a better estimator of the LG velocity in a O = 1.0 
than in a f2 = 0.3 Universe. As described earlier, in actual comparisons 
one uses the smoothed gravity of the LG as an estimator of its unsmoothet^ 
velocity. In Figure || we see that, for both models, the mean smoothed ve- 
locity of the LG candidates (dashed horizontal line) is smaller than their 
unsmoothed velocity, v = 600 km-s -1 . As stated in the previous Section, it 
means that part of the central velocity has been induced by mass fluctua- 
tions within the smoothing radius. Since (v R ) < v, though for the flat model 
(pR^ ~ (g R ) 5 {9 R ) should be smaller than v. Figure |I| show distributions of 
the values of g R . For f2 = 1, (g ) = 570 km • s _1 (dashed vertical line), so 
indeed (g R ) < v. Thus, the smoothed gravitational acceleration of the LG is 
in this case a slightly biased estimator of its unsmoothed velocity, and essen- 
tially all bias results from estimating unsmoothed velocity by its smoothed 
counterpart, and not from nonlinear gravity. The difference between (g R ) 



'Strictly speaking, smoothed over the small volume of the LG. 
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Figure 3: Modules of 500 km • s top-hat smoothed velocity and gravity 
for the LG candidates (scatter plots). Top panel: f2 = 1, bottom: Q = 0.3. 
Solid line represents the linear relation, dashed - mean values of g R and v R . 
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Figure 4: PDF of the smoothed gravity field (histogram). Solid vertical 
line marks the value of 600 km • s _1 , dashed - mean values of g R and shaded 
strips - la intervals around the mean values. 

and v results in a systematic error of the estimate of ft of 5.3 %. 

As for the model with = 0.3, the situation is different. The effects of 
nonlinear gravity act in the opposite direction than the effects of smoothing 
and turn out to dominate over them. The mean value of g R is here 642 
km • s _1 . In other words, in the $7 = 0.3 model the smoothed gravitational 
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acceleration of the LG overestimates its unsmoothed velocity. This results 
in a systematic error of the estimate of (3 of 6.5 %. 

To sum up, smoothed gravity of the LG underestimates its velocity in 
a high-fi Universe, while it overestimates the velocity in a low-fi Universe. 
Thus, there is no common bias, independent of the assumed model, that can 
be corrected for in the data analysis. Therefore, the results should be either 
expressed in a model-dependent way or they are a subject to this systematic 
error. Fortunately, this error is in any case relatively small, of a few per 
cent. 

Thus far, we have left aside the scatter in the velocity-gravity relation. 
This scatter results in a random error in the estimate of [3 from the LG 
velocity-gravity comparison. Unlike the systematic error, this error cannot 
be reduced in any model, because in practice one does not have at his 
disposal a whole set of the LG observers, but merely one observer (us). 
Thus, the NE put also an upper limit on the precision of this method of 
determining (3. The dispersion of the values of g , 5g R , is 62 and 79 km-s -1 , 
respectively for the models with = 1 and 0, = 0.3. This implies that the 
la random error in the estimate of (3 due to the NE is 12.3 % and 10.9 %, 
respectively. 

The summary of systematic and random errors is presented in Table EL 



n 


v/(g R ) - 1 


5g R /(g R ) 


0.3 


-6.5 % 


10.9 % 


1.0 


5.3 % 


12.3 % 



Table 1: Systematic and random errors in determining (3 from the LG dipole 
comparisons due to non-linear effects, for two cosmological models. 



5 Summary 

Using numerical simulations we have studied the relation between the ve- 
locity of the Local Group and its gravitational acceleration. This relation 
underlies both a test for the kinematic origin of the CMB dipole and a 
method of estimating (3. 

First, we have investigated the mean misalignment angle between the 
two vectors. For the case of Vt = 1, we have confirmed the result of Davis 
et al. (1991), that this angle is small. Actually, the value we have obtained 
(~ 7°) is even smaller than that of Davis et al. (~ 10°). We attribute a part 
of the difference to the different power spectrum we used, namely, that of the 
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PSCz survey. The rest of the difference is likely caused by the differences 
between the numerical codes (hydro vs. N-body). We have verified that 
the misalignment angle is fairly insensitive to £l, in agreement with our 
theoretical prediction. 

The observed misalignment angle is ~ 14°, close to the upper limit for 
the 95 % confidence interval for the simulated angle. Therefore, the NE are 
unlikely to be responsible for all of the observed misalignment. Other effects 
must play a role here as well, like, e.g., incomplete sky coverage and/or shot 
noise. 

Next, we have studied the relation between the amplitudes of the LG 
velocity and gravity vectors. This relation has a scatter, resulting in a 11- 
12 % random error of the estimate of (3. In the model with £1 = 1, the 
smoothed gravity of the LG turns out to be a biased low estimator of the 
LG velocity. On the contrary, in the model with Q = 0.3, the estimator is 
biased high. 

Using mock catalogs constructed for different cosmologies, S99 computed 
the ratio between the reconstructed gravity dipole at the observer's position 
and its true N-body velocity. In a 'numerical analysis', they subsequently 
used the average of all ratios as a multiplicative factor which relates the 
reconstructed gravity of the LG to the true LG velocity. Since biases due 
to the NE are distinct for different cosmologies (and this is also true for 
other effects, like, e.g., the effect of finite volume), an attempt to derive a 
model-independent estimate of (3 is illegitimate. Rather, the results should 
be expressed explicitly in a model-dependent way. Indeed, the estimates of 
(3 for models with £1 = 1 and with £l = 0.3, obtained by S99 using an alter- 
native, model-dependent maximum-likelihood analysis are clearly different. 
We conclude that the 'numerical analysis' of S99 is formally incorrect. 

While in principle important, in practice the discussed biases are small: 
of the order of a few per cent. With this accuracy, therefore, to the LG 
velocity-gravity comparisons the linear theory can be applied. Furthermore, 
the random error of the estimate of (3 due to the nonlinear effects is small 
compared to the total random error (~ 40%: S99). We thus conclude that 
in the LG velocity-gravity comparisons the NE are not the major concern. 
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